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A COMPARISON PRINCIPLE FOR THE POROUS MEDIUM 
EQUATION AND ITS CONSEQUENCES 


BENNY AVELIN AND TEEMU LUKKARI 


Abstract. We prove a comparison principle for the porous medium equation 
in more general open sets in R n+1 than space-time cylinders. We apply this 
result in two related contexts: we establish a connection between a potential 
theoretic notion of the obstacle problem and a notion based on a variational 
inequality. We also prove the basic properties of the PME capacity, in par¬ 
ticular that there exists a capacitary extremal which gives the capacity for 
compact sets. 


1. Introduction 

We study the porous medium equation (PME for short) 

f)v 

— - Au m = 0, (1.1) 

where m > 1. This equation is an important prototype of a nonlinear parabolic 
equation. The equation is degenerate, meaning that the modulus of ellipticity 
vanishes when the solution is zero. The name stems from modeling the flow of a 
gas in a porous medium: the continuity equation, Darcy’s law, and an equation of 
state for the gas lead to (1.1) for the density of the gas, after scaling out various 
physical constants. For more information about this equation, including numerous 
further references, we refer to the monographs [8] and [17]. 

The comparison principle is a fundamental tool in the theory of elliptic and 
parabolic equations. In particular, it can be used to define a class of supersolutions 
which is the counterpart for superharmonic functions in classical potential theory: 
we call a function a semicontinuous supersolution, if it satisfies the comparison 
principle with respect to continuous solutions. The definition is due to F. Riesz [15], 
and it makes the development of a nonlinear potential theory feasible. 

The comparison principle for parabolic equations is usually formulated for space- 
time cylinders, meaning sets of the form fix = 12 x (0 ,T). The boundary values are 
then taken over the parabolic boundary , where only the initial and lateral boundaries 
are taken into account. However, one often encounters situations where one would 
like to apply the comparison principle in sets which are not space-time cylinders. 
Thus our main objective is to establish a comparison principle for the PME in 
more general open sets in K™ +1 . Such a result is occasionally called the elliptic 
comparison principle, in reference to the fact that the time variable no longer has 
a special role. Moreover, the elliptic comparison principle can be used to develop 
the Perron method in general space-time domains, see [3, 12]. We also present two 
applications where such a comparison principle is indispensable. 
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For the heat equation, when m = 1, one may add constants to solutions. A com¬ 
parison principle for general open sets then follows from the space-time cylinder 
case by a straightforward exhaustion argument. For the PME, there is a compari¬ 
son principle over cylindrical domains, but adding constants is no longer possible. 
Our idea for circumventing this difficulty is to multiply one of the functions being 
compared by a constant close to one. The modified function is no longer a solu¬ 
tion, but it still satisfies the PME with an error-term on the right hand side. The 
error-term vanishes as the multiplicative constant tends to one. The comparison 
principle for the original functions then follows by the usual duality proof, modified 
to account for the error-term. Our argument yields a comparison principle for open 
sets of the form VLt n K, where K is a compact set. 

As the first application, we consider the obstacle problem. Roughly speaking, 
this amounts to finding a solution to a PDE subject to the constraint that the so¬ 
lution stays above a given function, the obstacle. Here we use a potential theoretic 
method for solving the problem: we define the solution to the obstacle problem to 
be the infimum of all supersolutions lying above the obstacle (reduite). For smooth 
enough obstacles the reduite is the smallest supersolution above the obstacle. The 
concept of reduite is standard in classical potential theory, and it has been uti¬ 
lized in a nonlinear parabolic context in [14]. Existence and uniqueness follow in 
a straightforward manner, at least for continuous obstacles. However, the relation 
between the smallest supersolution and the variational solutions to obstacle prob¬ 
lems constructed in [4] is not obvious. In this direction, we prove that the smallest 
supersolution is also a variational solution for sufficiently smooth obstacles. This 
follows from two facts. First, we prove that the smallest supersolution can always 
be approximated by variational solutions. Second, the notion of variational solution 
is stable with respect to the convergence of the obstacles in certain norms, see [4]. 
The converse of this, i.e. whether a variational solution agrees with the smallest 
supersolution, remains a very interesting open problem. 

The second application is a notion of parabolic capacity for the PME. This 
concept is defined via a measure data problem, as in [11] for the parabolic p- 
Laplacian. See also [18, 19] and the references therein for the capacity for the heat 
equation. We prove the basic properties of the capacity related to the PME, such as 
countable subadditivity and the existence of the capacitary extremal of a compact 
set. Our comparison principle plays a key role in the latter argument. 

The paper is organized as follows. In Section 2, we recall the necessary back¬ 
ground material, in particular various notions of supersolutions. Section 3 contains 
the proof of the comparison principle, and Section 4 is concerned with the obstacle 
problem. Finally, the basic properties of capacity are proved in Section 5. 

2. Weak supersolutions and semicontinuous supersolutions 

Let H be an open and bounded subset of M", and let 0 < t\ < t 2 <T. We use 
the notation fix = H x (0,T) and Ut lt t 2 = U x (ti,< 2 ), where U c fl is open. The 
parabolic boundary d p Ut lt t 2 of a space-time cylinder consists of the initial 

and lateral boundaries, i.e. 

d P U tl , ta = (Ux{t 1 })u(dUx[t 1 ,t2]). 

The notation Ut lt t a <i means that the closure U tl .t 2 is compact and U tl .t 2 c fl T . 

We use H 1 ( H) to denote the usual Sobolev space, the space of functions u in 
L 2 (S2) such that the weak gradient exists and also belongs to L 2 ( fl). The norm of 
is defined by 

IMlirHfi) = ll u lli 2 (n) + |V w llz, 2 (fi) ■ 
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The Sobolev space with zero boundary values, denoted by Hq (SI), is the completion 
of with respect to the norm of H 1 ( fl). 

The parabolic Sobolev space L 2 (0,T;7J 1 (f2)) consists of measurable functions 
u : Or [-o°, oo] such that x i-> u(x,t) belongs to for almost all t € (0,T), 


and 


f 

J fin 


|it| 2 + |Vu| 2 dxdt < oo . 


The definition of L 2 (0, T; Hq (f2)) is identical, apart from the requirement that 
x i-»- u{x,t) belongs to Hq( fl). We say that u belongs to L 2 )c (0, T ; Hf oc (U)) if 
u 6 L 2 (ti,t 2 ',H 1 (U)) for all U tl ,t 2 ^ 

Supersolutions to the porous medium equation are defined in the weak sense in 
the parabolic Sobolev space. 


Definition 2.1. A nonnegative function u : f It - 
equation 

^ - A u m = 0 
dt 

in n T , if u m e Lf oc (0, T; H^fl)) and 


is a weak supersolution of the 
( 2 . 1 ) 


[ 

J fin 


-u^- + S7u r ‘ 
at 


\Jpdxdt > 0, 


for all positive, smooth test functions ip compactly supported in f It- The definition 
of weak subsolutions is similar; the inequality is simply reversed. Weak solutions 
are defined as functions that are both super- and subsolutions. 


Weak solutions are locally Holder continuous, after a possible redefinition on a 
set of measure zero. See [7], [8], [9], [17], or [20]. 

We have also the following class of supersolutions. 

Definition 2.2. A function u : f lx -*■ [0, oo] is a semicontinuous supersolution, if 

(1) u is lower semicontinuous, 

(2) u is finite in a dense subset of CIt, and 

(3) the following parabolic comparison principle holds: Let d fl, and let 

h be a solution to (2.1) which is continuous in Ut 1} t 2 - Then, if h < u on 
dpU tl ,t 2 , h <u also in U tl ,t 2 - 

Note that a semicontinuous supersolution is defined in every point. Every weak 
supersolution is a semicontinuous supersolution provided that a proper pointwise 
representative is chosen. This is a consequence of the following lemma. 

Lemma 2.3 ([2]). Let u be a nonnegative weak supersolution to the porous medium 
equation in x (^ 1 ,^ 2 )- Then u has a lower semicontinuous representative. 


In the other direction, a bounded semicontinuous supersolution is also a weak 
supersolution, as shown in [13]. If unbounded functions are allowed, then the class 
of semicontinuous supersolutions is strictly larger, since the Barenblatt solution is 
a semicontinuous supersolution, but is not a weak supersolution, see [13]. 

Lemma 2.4 ([13]). Let u be a weak supersolution such that |u| < M < 00 . Then 


m|2 dxdt < 16M 2rra T 


IL t ^ u 

for all nonnegative functions e C™(Q). 


[jWl 

Jo. 


2 da; + 6M m+1 


f r/ 2 dx, 

Jn 


An application of the Riesz representation theorem shows that for each weak 
supersolution it, there exists a positive Radon measure p, u such that 

dp 


ff -u-^-+ \/u m ■ M <p dx dt = [ ipdp u 
JJn x at Jn x 
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for all smooth compactly supported functions p. This is the Riesz measure of 
u. The integrals on the left hand side do not depend on the particular point- 
wise representative of a supersolution. Thus a weak supersolution u and its lower 
semicontinuous regularization u have the same Riesz measures. 


Lemma 2.5. Ifu and v are weak supersolutions in Lloo, u,v - 0 on d p floo, u m ,v m € 
L 2 ( 0, oo; Hq(CI)), and p v < p u , then v < u a.e. in ^loo- 


Proof. Let p € C'^°(fl 0 o) be nonnegative. By subtracting the equations satisfied by 
u and v and using the assumption about the measures, we have 

/ -(it - v)<pt + V(u m - v m ) ■ Vpdxdt = I pdp u - / pdp v >0. 

S^oo ^ ^ rioo 


By a standard approximation argument using the fact that u m ,v m e L 2 (0, oo; Hq (fl)), 
we may also take the test functions p € C' 00 (fl 00 ) so that p - 0 on the lateral bound¬ 
ary <9f l x (0, oo). We apply the Green’s formula to get 



-(it - v)pt - (u m - v m )Apdxdt > 0. 


The fact that v < u follows from this inequality by repeating the standard duality 
proof for the comparison principle for the PME, see e.g. [7, Lemma 5], [8, Theorem 
1.1.1], or [17, Theorem 6.5]. □ 


Lemma 2.6. Let Ui, * = 1,2,... is a uniformly bounded sequence of weak superso¬ 
lutions in floo such that Ui -*■ u a.e. in Lloo- Then u is a weak supersolution in 
and 

lim / 4>dp Ui = / 4>dp u , 

*-*■«> ~/f2oo 


for every cj> e C“(fl<x>). 


Proof. Due to the uniform bound on the functions u,, it easily follows that 



d<t> 

- u m~ u 


A (j> dx dt > 0 . 


( 2 . 2 ) 


An application of Lemma 2.4 on each iq implies that Vu™ € L^ oc (fl 00 ). This, 
together with (2.2) yields that u is a weak supersolution. The claim about the 
measures follows from the computation 


lim f <j>dp u . = lim f -tij — + Vit™ • V</»dxdt 
»-» 00 -/Ooo Jflco ot 


[ 

f 

Jflc 


d(j> _ 

-u-£ + Vu” 
ot 

4>dpi u . 


V0dxdf 


□ 


We will frequently use the following characterization of the weak convergence of 
measures. See [10, Theorem 1, p. 54] for the proof. 

Theorem 2.7. Let p and pk, k = 1,2,3,..., be Radon measures on K™. Then the 
following statements are equivalent. 

(1) For all compactly supported smooth functions <f>, one has 

lim / (fdpk = / 4>dp. 

k ->oo Jr™ J 1” 

(2) For all compact sets K, one has 

lim sup pk{K) < p(K) ■ 

k->o o 
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(3) For all open sets U, one has 

/i([/) < liminf Hk(U) . 

k-* oo 


3. A COMPARISON PRINCIPLE 


The core of our arguments is a suitable form of the comparison principle, which 
we will prove in this section. We will work extensively with finite unions of space- 
time cylinders, so we begin by introducing some notation for such sets. For space- 
time cylinders Ut lt t? = U x (t 1; t 2 ), we denote the lateral boundary by 

S(U tl j 2 ) = dU x ■ 

For a cylinder the definition of the parabolic boundary is standard, but for finite 
unions of space time cylinders we will recall the definitions. The lateral boundary 
of a finite union of space time cylinders Uf 4 is then given by 

t l’ t 2 

S(u (u«S(^))x(uC^ 4 ). 

We also denote the tops of u Uh by 

r l’ r 2 

T(u^ 4 ):=(uCAx{4})x(u^ 4 ), 
and the bottoms similarly as 

(uC^x{ti})x(uC^ 4 ). 

Thus the parabolic boundary of Q = u[/\ ,,, is *S(Q) u B(Q) and the parabolic 

t l’ t 2 

boundary of backwards in time equations becomes S(Q ) u T(Q). 

We want to use the very weak (i.e. distributional) formulation of the porous 
medium equation, so we consider smooth test functions <f> e C°°(Q ) where Q = 
uUh i , such that <j) = 0 on S(Q). Note that the gradient of (j) does not necessarily 
vanish on S(Q). In the following we will always work with Cl a smooth domain. Let 
us now write the PME in terms of the above class of test functions. Assume at first 
that <f) has compact support in space. Then a standard approximation argument 
shows that we may write the definition of weak solutions as 

/ [ - ud>t + Vu m • Veil dxdt + / ud>dx- tub dx = 0 . 

JQ L J JT(Q) Jb(Q) 

After this, we may pass from compactly supported test functions to test functions 
vanishing on the sides *S(Q), since u and Vu m are in Lr(Q). Now, apply Green’s 
formula, which is justified by the usual trace theorem, to get 



- u m A(A] dx df 


+ f u m d n <j)dadt = 0 . 
JS(Q) 


(3.1) 


A similar argument can be carried out for weak supersolutions and subsolutions. In 
these cases, we get the appropriate inequalities in the final form. This formulation 
will be our starting point in the proof of the comparison principle. 


Theorem 3.1. Let K be a compact set in f It where f l is a smooth domain, let u 
be a nonnegative upper semicontinuous function which is a continuous weak super¬ 
solution in in fiy x K and satisfies u m e T 2 (0,T;Id 1 (fl)). Let v be a non-negative 
lower semicontinuous function which is a weak supersolution in f It x K and satisfies 
v m 6 L 2 (0,T■ H 1 (Cl)), v > 0 on K and u < v on K u 8 p CIt■ Then u <v in CLt . 
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Proof. We let e > 0, and denote 

D e = |(a :,t) 6 ^ ^ • 


The function 


■ v is upper semicontinuous, so that the set D t is closed in VLt- 


Moreover, the set D e does not intersect the set K , since u < v on K and inf k v > 0. 
Since K is compact, there is a positive distance between D e and K. Thus we can 
cover D e with a finite collection of space time cylinders not intersecting K. Denote 
the covering set by Df , and note that since D ei c D e2 for ei > e 2 we may choose 
the coverings for different values of e so that Df c Df. Let then Df = f It n Df. 
The set Df is still a finite union of space-time cylinders, and the function u is a 
weak solution in Df. 

Let u e = ff-. We want to compare u e with v in Df. To this end, note first that 
u e < v on dDf. Further, the function u e is a solution to 

(1 + e) 7 ™- 1 - 1 

(u £ ) t -AKr = /:=L_i__ A ^, 

interpreted in the sense of distributions. To see this, we compute 
[«e]t - A [uf\ m =K - «]t - A [u? - u m ] 


(3.2) 


(l + e ^ ((1 + e) m 1 ) 

1 


1 Au" 


\ 1 + e (1 + e) m ) 


Au m = 


(l + e) 


m— 1 


(l + e) 7 


-Au" 


We aim at adapting the proof of the comparison principle for the PME, see e.g. 
[7]. To proceed let 

D e.s = (O’A) 6 Df : t < s} , 

and take positive functions cj> e Cff(Df s ), and e C°°{Df s ) which vanishes on 
S{Df s ), d n if < 0 on S(Df s ) and so that if equals </> on T(Df s ). Denote also 
b = u e -v for brevity. Since b is negative and consequently also u™ - v m on dDf, 
we have from (3.1) and (3.2), since d n f < 0 on S(Df) 

[ fibdxdt = [ bchdx- f bibdx- [ bib t dxdt 

Jd f s Jr(Df a ) J B( Df s ) J Df s 

- f \ U T ~ v m ]Aijjdxdt + f [uf - v m ]d n ip dcr dt 
Jopf J S(Df 3 ) 1 

> f bfdx- f b’lptdxdt- f [uf - v m ]Aijjdxdt. 

Jt(d r.) Jd f. JDf. 


We can rewrite this as 

[ 

JT(D? S ) 

where 


J ~ ^ b(j)dx<J~ f bfift + aAijj) daidf + ^ fipdxdt, (3.3) 


0, 


if u e f v , 
if u e = v . 


Next we use a regularization to make the term 'ft + a.Af small in the above inequal- 

imooth fi 

<ci k <k , 


ity. To do this let Ofc, k - 1,2,..., be smooth functions in Df s such that 

1 
k 

and 

(a fc - a) 2 


Jd f 


■ dx dt -»• 0 as k 


CLk 


(3.4) 
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We replace the function ip in (3.3) by the solution ipk to the following boundary 
value problem 

Ut + ak Art = 0 , in D^ s , 

• u(x,s) = cp(x,s), on T(Df >s ), (3.5) 

u = 0, on S(D ^ s ), 

and get by Holder’s inequality 



a - cik)Aipk da; dt - 


L 


f'lpk dx d t 


< 


L 


b 2 (a Uk)2 dxdt 
ak 


1/2 r 


L 


au{Aipk) 2 dx dt 


1/2 


+ 



(3.6) 


To continue, we need to estimate the term on the right hand side of (3.6) contain¬ 
ing the quantity ak(Aipk) 2 independently of k. To do this we follow the calculations 
of [7]. We use the equation (3.5) for ipk and integrate by parts, first in time and 
then in space, which gives 


[a k Ail>k]/^4’k dxdt = - \ipk\t&ipk dxdt 

= / ipk[Aip k ] t dxdt- / 0A(pdx+ / ipkAipk dx 
JDf. JT(DF.) Jb(dF.) 


IT( D 1 S ) 

= f ipkA[ipk\t dx dt + f |V0| 2 dx - f \S7ipk\ 2 dx 
JDf s Jt(D^) ' Jb(DI„Y 

f d n ipk[ipk]t da dt + f Aipk\ipk]t dxdt + f |V</>| 2 dx. (3.7) 

JS(Df.) JDf. JT(Df.) 


Note now that the first term on the right hand side in (3.7) vanishes, since for 
almost every t < s we have [ipk\t = 0 on S[Df s ] due to the fact that ipk vanishes 
smoothly on the boundary. For the second term on the right hand side in (3.7), we 
use the first line in (3.7). This implies that 

r ak(Aipk) 2 dxdt < - f |V<(>| 2 dx. (3-8) 

JDf s 2 Jt(df s ) 

With the estimate (3.8) in hand, we see from (3.4) and the fact that b is bounded 
that 


Df. 


b 2 ( a ~ ak) 2 
ak 


dx dt ) (4 ak(Aipk) 2 dxdt j —> 0 as k-* 0. (3. 


9) 


To proceed we need to take care of the term involving / on the right hand side 
in (3.6). Recall that, as a distribution 

/ = K , } , --A u m . 

(1 + e) m 

Since the function ipk vanishes on the lateral boundary S(Df s ) of Z?f s , we have 


JDf. 


dx dt = 


(1 + 4 


-1 


(1 + e) 

(1 + e ) m_1 - 1 
(1 + e) m 


f Vu m • SJipk dxdt 
JDf. 


Df, 


|Vu m | 2 dx dt 


1/2 


Df 


\ 1/2 

|V'0fe| 2 dxdtJ . (3.10) 


By the assumption u m e L 2 (0,T;iJg(fl)), we see that the first integral is bounded 
independent of k and e. 
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Next we need estimate the L 2 -norm of |V"0fc| independently of k , which we do as 
in [17, p. 133]. Multiply the equation (3.5) for ip k by the test function 9 = AipkX(t), 
where \(0) = 1/2 and x(s) = 1; thus Xt “ y Next we integrate by parts, first in 
space and then in time, and get 


0= / [ip k ] t Ai/j k xdxdt + / a k (Ai/; k ) 2 xdxdt 

Jdf s 

= - / V(V , fc)« ' W’fcXdxdt + / a k (Ai/j k ) 2 xdxdt 

JDf s 

= ~\ f [\S 7 -tpk\%xdxdt + f a k (Aip k ) 2 xdxdt ( 3 . 11 ) 

=5 d.Td( - i J t(du wux 

+ ^ f, p J v ^/0| 2 x+ [ „ a k (Aip k ) 2 xdxdt, 

2 Jz3(D e F J Jd f s 

using (3.5) and (3.11) we get 


- [ (]VV , fc| 2 )decd t+ f a k (Ai/; k ) 2 xdxdt 

s JD? a Jd£, 

\Vip k \ 2 Xdx- [ \Vip k \ 2 xdx\<C [ IV</>| 2 da;. (3.12) 

Jb(d£,) ] Jt(df s ) 


<C 

' >r(ni,) 

Combining (3.6),(3.9),(3.10) and (3.12), we have so far established 

, (1 + e) m_1 - 1 


fr 


(u e - v)d)dx < C- , 

Ji (1 + ey 


fr 


T(Df.) 


]V</>r dx. 


(3.13) 


Before letting e -»• 0, we still need to check that 



dx <C, 


for some constant C not depending on e > 0. We are free to assume that </> e 
C™(D t0}S ) for some eo- Then, since -D £0jS c Df s for e < eo, we have T(Df' s )nD e0tS c 
T(D e0}S ) which proves the desired bound. Thus, letting e -»■ 0 in (3.13), we get that 



n[R n x{s}] 


(u - v)cj)dx < 0. 


Since this holds for any positive </>, we obtain that u < v a.e. in fly n [K" x {s}] for 
any s, and then also in Qt- □ 


The crucial point in the proof above is that we can approximate the set 
D e = j(x,f) 6 fl T ■ ^ > 'cj , 

by finite unions of space time boxes while staying inside the set where u is a weak 
solution. Thus we can also deduce the following Theorem: 

Theorem 3.2. Let E be an open set in R” +1 , let u be a non-negative continuous 
weak solution in E such that 

[ [\u m \ 2 + \Vu m \ 2 ]dxdt < OO . 

Je 

Let v be a non-negative lower semicontinuous weak supersolution such that 

f [\v rn \ 2 + |Vv m | 2 ] dxdt < oo, 

JE 
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v > 0 on dE and u < v on dE. Then u < v in E. Furthermore if a connected 
component of dE is the boundary of a finite union of space-time cylinders then we 
can remove the assumption v > 0 on that component. 

4. The obstacle problem 

In this section, we construct solutions to the obstacle problem by a potential 
theoretic method. More specifically, we call a function u a solution to the obstacle 
problem if it is the smallest supersolution above the given obstacle function ip. 

Existence and uniqueness are fairly easily established for this notion of solution to 
the obstacle problem. However, the relationship between the variational solutions 
studied in [4] and the smallest supersolution is not immediately clear. In this 
direction, we apply the comparison principle established earlier to prove that the 
smallest supersolution is also a variational solution, provided that the obstacle 
is sufficiently regular. This is a consequence of two facts: first, we prove that 
the smallest supersolution is a point-wise limit of variational solutions. Second, 
variational solutions are stable with respect to convergence of the obstacles in a 
suitable norm. 

We expect that the converse is also true, i.e. that a variational solution is the 
smallest supersolution . However, our version of the comparison principle in general 
domains is not strong enough to prove this. 

First we describe the notion of smallest supersolution in more detail. 

Definition 4.1. Let ip be a positive, bounded measurable function in and 
denote 

U ^ = {u is a semicontinuous supersolution in : v >ip in floa} . 

We define the reduite (or reduced function) of ip as 

I?,/, = inf{v : v e U^} . 

For a measurable set E, we abbreviate Re = 1? XE . We denote by R^, (lower semi¬ 
continuous) ess lim inf-regularization of R^. The function R.^, is usually called the 
balayage of ip. 

The terms reduite and balayage come from classical potential theory. The notion 
is due to Poincare. We will need the following basic theorem, for which the proof 
is standard, but we reproduce it here for the reader’s convenience. 

Theorem 4.2. The balayage R^ is a semicontinuous supersolution in Ht- 

Proof. Pick a space-time cylinder ^ and a weak solution u which is con¬ 

tinuous in Ut lt t 2 with u < R^ on d p Ut lt t 2 • Then also u < v on d p Ut 1 ,t 2 for v e U 
and by comparison the same holds in Ut 1 ,t. 2 - We take the infimum over v to get 
that u < R^. Since u = u by the continuity of u, we conclude that u < R^. □ 

Note that in general, R^ might not be lower semicontinuous, and R^, might 
not be above ip in every point. However, for continuous ip it holds that R^ > ip 
everywhere. This together with Theorem 4.2 implies that R^, is the unique smallest 
semicontinuous supersolution above the obstacle ip. By the smallest supersolution, 
we mean a function u e U ^ with the property that 

u<v for all veU^j. (4.1) 

A semicontinuous supersolution with the property (4.1) is unique, if it exists; in¬ 
deed, if there are two functions U\,U 2 e lAf satisfying (4.1), then two applications 
of (4.1) give the inequalities u i < U 2 and U 2 < iti, so that u\ = U 2 - 
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The next aim is to relate the smallest supersolution to the variational solutions 
to the obstacle problem constructed in [4], We first recall some facts from [4], 

We consider nonnegative obstacle functions ip defined on f It, with compact 
support and satisfying 

1> m €L 2 (0,T-,Hl(n)), d t (ip m )eL B &(Sl T ). (4.2) 

The class of admissible functions for the obstacle problem is defined by 

Ify(fir) := {v ■■ n T [0,oo] : v m e L 2 (0,T;,v >ip a.e. on U T }. 

Note that ip 6 and therefore + 0. 

With the above classes, we can state the definition of a strong solution to the 
obstacle problem. 

Definition 4.3. A nonnegative function u e is a strong solution to the 

obstacle problem for the porous medium equation if dtu e L 2 (0,T;iJ _1 (fl)) and 

f T (d t u,a(v m -u m ))dt+ f aVu m ■ V(v m - u m ) dz >0, 

J 0 J fix 

holds for all comparison maps v e K^>(Qt) and every Lipschitz continuous cut-off 
function a : [0,T] -*• [0,oo] with a{T) - 0. 

The cutoff function a is needed for making this definition consistent with the 
definition of weak solutions to the obstacle problem, which we will recall later. 

For the existence of strong solutions, we still need the assumption 

^:=d t V>-AV> m eL°°(Q T ). (4.3) 

The following result can be extracted from [4, Theorem 2.6]. 

Theorem 4.4. Let Q be a bounded open subset of R" with a smooth boundary. 
Assume that the obstacle ip satisfies the regularity conditions (4.2) and (4.3). Then 
there exists a strong solution u to the obstacle problem for the PME in the sense of 
Definition f.3 satisfying u m e L 2 (0, T ; Hq(Q)) and u( ■ ,0) = 0. 

The function u is also locally Holder continuous, and satisfies u>ip everywhere 
in f It- Further, u is a weak supersolution to the porous medium equation in Qt, 
and a weak solution in the open set {z e 0,? '■ u(z) > ip(z)}. 

We now wish to show that u in Theorem 4.4 is a weak solution in the larger set 
[fir x supp(^)] u {z € fl T : u{z) > ip(z)} . 

With this in mind we recall the following form of a partition of unity. 

Lemma 4.5 (Partition of Unity). Let U\, U 2 , ■ ■ ■, U n be open sets, and let K be a 
compact set such that K c U\ u U 2 u ••• u U n . Then there exist functions rji 6 Cff(Ui) 
such that 

n 

Yj Vi = 1 0Tl K • 

i=l 

Proof. For a version where the functions r)i are continuous, see [16, Theorem 2.13, 
p. 40]. The fact that one may also choose smooth functions follows easily from the 
continuous version by applying a suitable mollification. □ 

Lemma 4.6. The strong solution to the obstacle problem given by Theorem 4.4 is 
also a weak solution to the PME in the set Qt x supp(^>). 

Proof. Let S > 0 be a number, and let rjs ■ R [0,1] be a Lipschitz function with 
r/(s) = 0 for s < -5, ry(s) = 1 for s > 0, and |r/(s)| < 1/S. The solution u is constructed 
in [4] as the uniform limit as S -*■ 0 of solutions to 

d t us - A uj 1 = r, 5 (ip m - u?)(d t iP - A if m ) + . 
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The claim now follows from the fact that (dtip - A ip m )+ - 0 in VLt n supp(i/)). □ 

Theorem 4.7. Let ip be a nonnegative, compactly supported function satisfying the 
regularity assumptions (4.2) and (4.3), let u be the strong solution to the obstacle 
problem given by Theorem 4-4 with obstacle ip, and denote K = supp(i/') n {u = ip}. 
Then u is a weak solution in floo x K. 


Proof. We first show that the function it is a weak solution in \ K. Denote 
U\ = floo n supp(?/;) and U 2 = {u > 4’}. These sets are open, and 

(loo \ K = Ui u U 2 ■ 


Further, it is a weak solution in U\ and in Ui- The claim concerning the set U\ is 
Lemma 4.6, and the claim about U2 is a part of Theorem 4.4. To show that u is a 
solution also in U\ u C/ 2 , let e u[/ 2 ). An application of Lemma 4.5 shows 

that there are functions rji € Cpf (Ui), i = 1,2, such that 771 + 772 = 1 on the support 
of ip. By applying the fact that u is a weak solution in Ui and in U2, we get 

f -udtip + Vu m ■ \7pdxdt = V f -udt(<pi]i) + Vu m ■ Vdxdt = 0. 


Since this holds for any test function ip, u is a weak solution in U\ u f7 2 - 


□ 


We are now ready to proceed with the approximation result. 

Theorem 4.8. Let ip be continuous and compactly supported in Qt- Then the 
smallest supersolution R^> is an increasing limit of strong solutions Wj to the ob¬ 
stacle problem with smooth compactly supported obstacles (pj increasing to ip. 

Proof. Let Uj = {4’ > 1 /j} for j = 1,2,..., Kj = Uj. First note that Kj and Kj+i 
have a positive distance between them. Thus if we let hj = (\/ip - we see 

that hj+\ - hj is strictly positive in Kj. By a mollification argument it can easily 
be seen that for each j = 1,... there exists a function fj e Cpf (A' ?+1 ) such that 

hj < fj < h j+ i. 

Taking cpj = /?, j = 1,..., we immediately see that <p™ e Cq (Kj+i), thus it satisfies 
(4.2) and (4.3). Moreover by construction we get 

<pi < cp2 < ■ ■ ■ < ipi 4>j ip as j 00 • 

Let Wj be the strong solutions to the ^-obstacle problems. Since > ip by the 
continuity of ip, we have Wj < R^, on K = d({wj = <pj} n supp^j). Also note that 
K c Uj+i , whence R^, > —j- > 0 on K. This allows us to use the comparison principle 
of Theorem 3.1 together with Theorem 4.7 to get that w 3 < R^,. A similar argument 
shows that Wj < Wj+ 1 . Thus w = Hindoo Wj is a semicontinuous supersolution as an 
increasing limit of continuous supersolutions, and w < R^. To finish the proof, we 
have that w > ip everywhere in £It, whence R^, < w. Thus 

w < Rj, < Rj, < w , 

and the proof is complete. □ 


The final step is to combine the approximation result with a stability result for 
variational solutions to conclude that the smallest supersolution is also a varia¬ 
tional solution. We recall some more facts from [4], in particular the notion of a 
weak variational solution, for which stability with respect to the obstacles can be 
established. 

For the notion of weak solutions, we use the class of admissible comparison 
functions 

= {v 6 K^Qt) : d t (v m ) € L^(n T )} . 
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We need to make sense of the time term in the variational inequality when we do 
not know that d t u belongs to the dual of the parabolic Sobolev space. We do this 
as in [1] and [4]. We recall the notation 

L ->[4 1 

J L L i 


« d t u,a V (v m -u m ))) U0 = 


+q(o) X’'[ 


TO + 1 
1 


-u m+1 - uv r - 


- audtv r ‘ 


dx df 


to + 1 


0 ) 


da; 


where uq e L m+1 (Q) is a function giving the initial values of the solution, and a 
is a nonnegative Lipschitz continuous cutoff function depending only on the time 
variable with a(T) = 0. The role of the function a is to eliminate the final time 
term, as we do not know in general whether u is continuous in time. Observe that 
if d t u e T 2 (0,T; ff _1 (f2)), we have 

[ T (d t u , a V (v m - u m ))dt = ((dtu, a V (v m - u m ))) UQ . 

J o 

This follows formally from integration by parts, and the rigorous justification is 
given in [4, Lemma 3.2]. This makes the following definition consistent with the 
definition of strong solutions in the previous section, i.e. strong solutions are also 
weak solutions. 


Definition 4.9. A nonnegative function u e is a weak solution to the 

obstacle problem for the porous medium equation if the inequality 

(( d t u,ar,(v m -u m ))) U0 + f aVu m -u m ))&z> 0 

holds true for all comparison maps v e and every nonnegative, Lipschitz 

continuous cut-off function depending only on the time variable with a(T) = 0. 

Theorem 4.10. Let LI be a bounded open subset of W 1 with a smooth boundary. 
Assume that the obstacle if satisfies the regularity condition (4.2). Then there 
exists a weak solution u to the obstacle problem for the porous medium equation 
in the sense of Definition f.9 satisfying u m e L 2 (0,T;Ug(fl)). Again, u is a weak 
supersolution to the porous medium equation in LIt- 

The following theorem may be extracted from the proof of Theorem 2.7 in [4]. 

Theorem 4.11. Let if, be a sequence of obstacles satisfying (4.2) with compact 
support in Qt such that 

V>r^V’ m inL 2 (Q,T-Hl(Ll)), and d t (ifT) - d t (if m ) in L^ (ft T ) , 

furthermore let Ui be the respective variational weak solutions to the obstacle problem 
with obstacle if,, see Theorem f.10. 

Then there is a function u e L°°(0,T; L m+1 (Ll)) with u m e L 2 (0, T; Hq(LI)) and 
u(-,0) =0 such that, up to subsequences, 

Ui u a.e., U™ -*■ u m in L 2 (LIt), and Vu™ -*■ Vu m weakly in L 2 (LIt) ■ 

Furthermore, u is a variational weak solution to the obstacle problem with obstacle 
if and initial values zero. 

Since strong variational solutions are also weak variational solutions, we get the 
following theorem as an immediate consequence of Theorems 4.8 and 4.11. 

Theorem 4.12. Let if be a continuous function with compact support in f It satis¬ 
fying the regularity assumptions (4.2). Then the smallest supersolution R^, is also 
a variational weak solution. 
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A general converse for Theorem 4.12 remains open. We record the following 
partial result for use in Section 5. 

Theorem 4.13. Let ip be a smooth obstacle with ip > 0 in Qt- Then any variational 
strong solution u to the obstacle problem coming from Theorem 4-4 satisfies u = R^. 

Proof. Since ip > 0, any strong solution is strictly positive inside fir- Thus, given a 
semicontinuous supersolution v e U^, we may apply Theorem 3.1 on the set {u > ip} 
to conclude that u <v. Since u e U^,, we get u = R ,/,. □ 


5. Parabolic capacity for the porous medium equation 

In this section, we define the parabolic capacity for the porous medium equation 
and establish its basic properties. 


Definition 5.1. The PME capacity of an arbitrary subset E of Qoo is 

cap(E) = sup{/.i(n oo ) : 0 < u,^ < 1, supp(/i) c E} , 

where p is a positive Radon measure, and u^ is a weak supersolution with = 0 
on d v Qoo, and a weak solution to the measure data problem 

Our next result is that there exists a capacitary extremal for the PME capacity of 
a compact set K, i.e. a semicontinuous supersolution u such that cap(AT) = p u (K). 
We need the following two lemmas. 


Lemma 5.2. Let ip be a smooth, positive compactly supported function, and set 
ip e = (ip m + £ m y/ m an d v e = R^. 


Then the limit function 


v = lim v F 


is a continuous weak supersolution and a weak variational solution to the obstacle 
problem with obstacle ip in Qt- Further, v is a weak solution in the open set {v > ip}. 


Proof. The existence of the point-wise limit as e -> 0 follows from the fact that 
Ri/j ei < R-ip e2 if Si < £ 2 - The limit v is an upper semicontinuous weak supersolution 
as a decreasing limit of continuous weak supersolutions, and v > ip since v E > ip E . 
By Theorem 4.13, we may take v E to be a strong variational solution to the obstacle 
problem. Hence v is a weak variational solution to the obstacle problem by [4] . The 
continuity follows from [5]. 

Since v E is a variational strong solution to the obstacle problem, it is weak 
solution in the set {v E > ip E }. If K is now a compact set contained in {v > ip}, we 
have that K is also contained in {u e > ip E } for all sufficiently small e, since 

v E -ip E > inf(u -ip) - e 

by the inequalities v E > v, -ip E = ~{ip m + £ m ) 1 / m > -ip-e, and the fact that inf/^(u- 
ip) > 0. Thus 

I -vd t ip + Vu m • Vipdxdt = lim / -v E d t ip + SJv™ ■ Vtpdxdt = 0 
Jn T e->oJn T 


for all smooth test functions ip with support in K. Since K was arbitrary, v is a 
weak solution. □ 
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The next lemma is the key step in constructing the capacitary extremal. For the 
proof, we record the following estimate. Let u be a positive weak supersolution in 
floo, u vanishing on the lateral boundary of floo- Suppose in addition that there 
exists a time to > 0 so that u is a weak solution in 12 t0jOo . Then 

u(x, t) < c(t - (5-1) 

for all t > to with a constant depending only on n, to, and the diameter of 17. This 
is the so-called universal estimate. See Proposition 5.17 in [17] for the proof. 

Lemma 5.3. Let K be a compact subset of Assume that u and v are lower 
semicontinuous weak supersolutions in lloo and that u is continuous in floo, and a 
weak solution after a time T such that I\ d f It- Moreover, assume that u > 1 in K, 
u = 0 on dpLloo, 0 < v < 1 in and v = 0 on <9 p 17oo- Then 

jl U (Ll oo) • 


Proof. Let 0 < ipi e Co°(^°°) be an increasing sequence of smooth obstacles con¬ 
verging to v such that ipi < v and ipi < ip%+\. Denote the perturbed obstacles 
ipi = (ip™ + e™) 1 !™, and the corresponding solutions to the obstacle problem by 
v\. Further, let Vi = lim^o^i be the weak supersolutions in 12oo constructed in 
Lemma 5.2. We argue as in [13, proof of Theorem 3.2, p. 148-149] to see that 
Vi < Vi+i < v and Vi -*■ v as i -*■ oo. Finally, we have that 

sup v\ = sup ipi < 1 + e 

f2oo 


and 


sup Vi = sup ipi < 1. 

i~2oo ^oo 

We define the supersolution 

w\ = min(u-, u) . 

By lower semicontinuity, the set {u > 1} is open, and K is compactly contained 
in it. This allows us to construct a compact set K' and an open set U such that 
K c U c K' c {u > 1}. If e is small enough we know that 1 + e < u in K', so that 
w\ - v\ in K'. Hence for such e we have for <p' = 1 on U and <p' = 0 outside K' that 


(U) < (p' dp v e = f <t>' dp w * ^ fJ-w f ( K'). 


(5.2) 


Next note that since u = 0 on d p Ll and u < e for sufficiently large times by (5.1), 
there is a compact set K" d K' in 12 c*, such that in \ K" we have wf = u since 
v\ > e. . Hence we obtain for <p" e Cff(Ll <*,) such that <p" = 1 on I\" that 


[ (P"dp w c = f 
Jfleo Jfl 

-[ 
J a 


M" 


U)-—- + V(wl) m • V 4>" dx dt 
at 


df_ 
dt 

f <P"dpi u . 


u~— — \-S7u™ • V(p dxdt 
at 


Thus we obtain the estimate 


Pw‘(K') < / <p" dfi W € = / cp" dp u < /^„(12oo) • 

D w )j OO ^ OO 


We combine (5.2) and (5.3) to get the inequality 


Hv‘(U) < Pu(^-oo) ■ 


(5.3) 
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By construction vf -*■ Vi point-wise, thus from Lemma 2.6 we get that p v « -»■ p Vi 
weakly. By the standard properties of weak convergence of measures, see Theorem 
2.7, we get that 

p Vi (U) < liminf/x„e(?7) < p u (Cloo) ■ 

£-*•0 ’ 

The sequence ( t>i ) is increasing, and converges point-wise to the original supersolu¬ 
tion v. Again from Lemma 2.6 we get the weak convergence of the corresponding 
measures. Another application of Theorem 2.7 now shows that 

Pv (A) < p v (U) ^ lim inf Pv i ( U'} < p u (floo), 

i—> oo 

and the proof is complete. □ 

A consequence of Theorem 3.1, is that in the special case that we have a de¬ 
creasing sequence of smooth obstacles converging to a characteristic function of a 
compact set, the obstacle problem is stable. If we had a full elliptic comparison 
principle, this lemma would hold for a decreasing sequence of smooth obstacles 
converging to an upper semi-continuous obstacle. 

Lemma 5.4. Let K c be a compact set. Let Ei d Lloo, i = 1,... be a shrinking 
sequence of open sets such that Ei +1 d Ei 

OO 

n Ei = K. 

i =1 

Assume that the non-negative functions ifi : -» M are supported in Ei, satisfy 

(4.2) and (4.3), and ifi > Xk, i = 1,2,..., is a decreasing sequence such that ipi -> \k 
point-wise in 0 as i -* oo. Then R x j, i -*• Rk point-wise in fioo and -*• pr k 
weakly as i -»• oo. 

Proof. By Theorem 4.7, the functions R^ )i are continuous. Thus an application 
of Lemma 2.6 shows that u - Hindoo R^. is an upper semicontinuous weak super¬ 
solution, and the respective measures also converge weakly. Further, 

u > Rk , 

since R^ Ji > Rk for each i. 

The lemma now follows if we prove the opposite inequality. To this end, note first 
and that from Theorem 4.7, R^ is a weak solution in {R^ i > (floo \supp(i/’i)), 
so that the support of the measure pr^. is contained in supp^j) c Ei. These sets 
shrink to K, and the measures pr^. converge weakly to p u . Thus supp(/z u ) c K, 
which implies that u is a weak solution in Id,*, \ K. 

If now v > Xk is an arbitrary semicontinuous supersolution with v = 0 on dpfloo, 
it follows from Theorem 3.1 that u < v. We take the infimum over v to get that 

u < Rk , 

and the proof is complete. □ 

A consequence of the stability Lemma 5.4 is that we have stability of the balayage 
with respect to decreasing sequences of compact sets. 

Lemma 5.5. Let Ki c i = 1,2,..., be a decreasing sequence of compact sets 
and denote K = n^Ki. Then Rk, is a decreasing sequence converging to Rk, 
moreover p^ converges to Pk k , weakly as i -»• oo. 

Proof. Let us construct £) = {d((x,t); Ki) < c/i}, i = 1,..., for a small constant 
c < 1 such that E\ c then the sequence Ei satisfies the requirements of Lemma 
5.4. 
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Let us now construct smooth functions e Cff(Ei) such that ifi = 1 on Aj, 
then let ifi = [if] 2 , and we have that R ki > Rk, by construction. As in the proof of 
Theorem 4.8, the sequence ifi will satisfy (4.2) and (4.3). It is now clear that the 
sequence ipi satisfies all requirements of Lemma 5.4 and thus we get that R ki -*■ Rk 
and consequently also Rk, -*■ Rk, furthermore using Lemma 2.6 we see that the 
measures Pr k . converge weakly to g.R K as i ->■ oo. □ 

Theorem 5.6. Let K be a compact subset offloo- Then 

cap (K) = h Ak (K). 

Proof. Since Rk is a semicontinuous supersolution such that 0 < Rk < 1, it follows 
immediately from the definition of the PME capacity that 

Tr k ( K ) ^cap(A), 

since Rk is a solution outside K. 

To prove the opposite inequality, let first K' d fbe a compact set such that 
K is K'. To be able to use Lemma 5.4 we will let Ei c K', i = 1,... be a shrinking 
sequence of open sets such that 


n Ei=K. 


2=1 


Let ifi e C™ (Ei), i = 1,..., be a decreasing sequence of smooth functions converging 
to xk point-wise in fas i -*■ oo, and such that 


ifi = 


1 + — 

2 * 


on K. 


Consider now the functions ifi = [V’,:] 2 , then if™ e C^Ei), and it is a decreasing 
sequence of functions converging to \k point-wise in fas i -*■ oo, such that 




on K , 


moreover ifi satisfies (4.2) and (4.3) for all rn > 1. Denote by Ui the corresponding 
solutions to the obstacle problems with obstacle ifi. Let now v be a weak superso¬ 
lution in floo such that 0 < v < 1 and r = 0on dpfloo- Then it follows from Lemma 
5.3 that 

I^V ( H ^ t^Ui (floo ) — P'Ui (R ) ■ 

We use Lemma 5.4 to see that p Ui -* weakly. The claim now follows from the 
above estimate, since 

lim sup fj, Ui (K') < n k]< (K') = ^ (K) 


by Theorem 2.7. 


□ 


We have now developed all the technical tools needed to establish the basic 
properties of the PME capacity, including that it is a regular, subadditive capacity. 

Theorem 5.7. The PME capacity has the following properties. 

(1) Countable subadditivity: In other words if Ei, i = 1,2,..., be arbitrary 
subsets office and E = u °fiEi, one has 

oo 

cap(E) < ^cap(A). 
i=1 
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(2) Stability with respect to increasing sequences of sets: Let E it i - 1,2,, 
be arbitrary subsets of floo with the property E\ c E 2 c •••. and denote 
E = u™ 1 Ei. Then 

lim cap (Ei) = cap(A). 

i-* 00 

(3) Stability with respect to decreasing sequences of compact sets: Let Ki c 
Qoe, i- 1,2,..., be a decreasing sequence of compact sets and denote K = 
n“i Ki. Then 

lim cap (Ki) = cap(iv). 

i—> 00 

(4) Let U d floo be an open set. Then 

cap (U) = p Ru (fl OO ) • 

Proof. From the methods developed in [11] we see that (1) and (2) follow from 
Lemma 2.5. Property (3) is a consequence of Theorem 5.6 and Lemma 5.5. Property 
(4) follows from (2), Theorem 5.6, and Lemma 2.6 as in [11, Lemma 5.9]. □ 

In conclusion we have established more than enough to say that Borel sets are 
Choquet capacitable: 

Theorem 5.8. The PME capacity is Choquet capacitable (inner regular). This 
means that for all Borel sets E c it holds that 

cap (E) = sup{cap(A'): K c E,K compact } . 

Proof. Since the capacity is monotone, stable with respect to increasing sequences of 
sets (Theorem 5.7, (2)) and stable with respect to decreasing sequences of compact 
sets (Theorem 5.7, (3)), it is a regular capacity and hence the claim follows from 
Choquet’s capacitability theorem [6, Theorem 9.3, p. 155]. □ 
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